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We have studied the nature of the surface charge distribution in CeTes- This is a simple, cleav- 
able, layered material with a robust one-dimensional incommensurate charge density wave (CDW). 
Scanning tunneling microscopy (STM) has been applied on the exposed surface of a cleaved single 
crystal. At 77 K, the STM images show both the atomic lattice of surface Te atoms arranged in a 
square net and the CDW modulations oriented at 45° with respect to the Te net. Fourier transform 
of the STM data shows Te square lattice peaks, and peaks related to the CDW oriented at 45° 
to the lattice peaks. In addition, clear peaks are present, consistent with subsurface structure and 
wave vector mixing effects. These data are supported by electronic structure calculations, which 
show that the subsurface signal most likely arises from a lattice of Ce atoms situated 2.53 A below 
the surface Te net. 

PACS numbers: 61.44.Fw, 68.37.Ef, 71.15.Mb, 71.45.Lr, 72.15.-v, 73.20.-r, 73.20.At 



peaks in the Fourier transform of the STM images were 
identified as satellite structures, consistent with the dis- 
commensuration picture. The data presented here are of 
higher quality and are compared to a compelling wave- 
vector-mixing analysis, originally suggested by Fisher et 
al. [10, 11], which does not involve discommensurations. 

The undistorted crystal structure of CeTes is shown in 
Fig. [H It is of NdTes p^, H^] type, weakly orthorhombic 
and described within the space group Cmcm. It is a 
layered structure that consists of two building blocks: 
double layers of [Te]~ square-nets, and puckered ionic 
[Ce2^Te2~]^^ double layers that are placed between the 
nets. The three-dimensional structure is composed of 
slabs of these structural motifs. The atoms within slabs 
are covalently bonded, while bonds between the slabs are 
weak van der Waals type, allowing the crystals to cleave 
easily between the Te layers. Hence the exposed surface 
is the Te-net, ideal for STM studies of the CDW. 

The existence of a unidirectional CDW in tritellurides 
was first reported in a transmission electron microscopy 
study of a series of RETea crystals by DiMasi and col- 
laborators [l[. They identified superlattice reflections in 
the electron diffraction pattern corresponding to a single 
incommensurate modulation wave vector with a magni- 
tude of qcDW ^ 2/7 X 27r/c, where c=V^ao and ao is the 
Te-Te separation of 3.1 A. This indicated the presence of 
incommensurate distortions in the Te layer planes. The 
incommensurate superstructure was solved within space 
group C2cm(007)000 for RE = Ce, Pr, and Nd by Malli- 
akas and collaborators using single crystal x-ray diffrac- 
tion |3|]- This study quantified distortions in the Te net, 
and revealed that a distribution of planar Te-Te distances 
exists, with a minimum value of ~ 2.95 A and a maxi- 
mum value of ~ 3.24 A. 

Due to the large lattice constant along the b-direction 

25 A), as shown in Fig. [TJ the Brillouin zone of 

* Electronic address: ICorresponding author. E-mail: tessmer@pa.msu.ed5] ETe3 is COmpressed (plate-like) and slightly orthorhom- 



I. INTRODUCTION 



The metallic compound CeTes belongs to a family of 
layered RETea materials, where RE is a rare-earth el- 
ement, which have received much attention as a model 
system to study incommensurate charge density waves 
(CDWs) [1, 2, 3]. This class of materials features two- 
dimensional square- net motifs composed of Te atoms. 
Such square-net arrangements have been considered the- 
oretically performing electronic band structure calcula- 
tions and were found to be susceptible to CDW formation 
driven by Fermi surface nesting [4] . On the experimental 
side, measurements have shown rather weak coupling be- 
tween the layers and large energy gaps as high as 400 meV 
for CeTes [5|,i,[7]. The CDW is well formed in CeTes at 
room temperature, and no transition to a non-CDW state 
has been observed for temperatures as high as 500 K [8] . 
Although it is well established that the CDW forms in 
the Te net, the exact nature of the CDW in the RETes 
family has not been resolved to date. For example, there 
is an ongoing debate regarding whether the CDW is uni- 
formly incommensurate or locally commensurate within 
domains, with phase slips, i.e. discommensurations, oc- 
curring at the domain walls [1, [l^ . 

In this paper we present scanning tunneling microscopy 
(STM) measurements obtained at a temperature of 77 K 
and theoretical calculations of the CeTes surface, with 
the main focus on understanding the surprisingly-large 
subsurface contribution to the tunneling signal. This 
study can be compared to room temperature STM data 
presented in Ref. [91]. Ref. also included Pair Dis- 
tribution Function analysis of x-ray diffraction data that 
showed clear evidence for discommensurations; moreover, 
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FIG. 1: The average (undistorted) crystal structure of CeTes 
consisting of corrugated CeTe slabs, and Te layers, where Te 
atoms are separated by 3.1 A in a square- net. The side view 
shows two complete unit cells, with the unit repeated along 
the b axis. The top view shows the surface Te net exposed 
upon cleaving the crystal. The first subsurface layers of Ce 
and Te atoms are indicated. Taking into account atomic re- 
laxation, with respect to the b-direction (perpendicular to the 
layers), the subsurface Ce lattice is 2.53 A below the surface 
Te net; the subsurface Te lattice is 3.44 A below the Te net. 



FIG. 2: A real-space STM image of the Te net obtained at 
77 K, showing both Te atoms and CDW modulations oriented 
at 45° to the net. The image is the average of four images that 
were obtained consecutively, in constant-current mode, at a 
at a sample bias of 100 mV and tunneling current of 0.6 nA; 
no image processing or filtering of the data was performed. 
Here we show the largest area without substantial contam- 
ination, although some is present, as seen in the upper left 
corner. The lines on the lower right indicate locations of high 
charge density due to the CDW, while the arrow marks the 
CDW wave- vector direction. The expanded image on the up- 
per right includes a square grid as a guide to the eye; a Te 
atom is located at each intersection. 



bic [1^. Strong anisotropy has been observed in their 
transport properties [Jj, |6|] , which reflects the nearly two- 
dimensional nature of the system originating from the 
weak hybridization between the Te layers and the RETe 
slabs. Electronic structure for tellurium planes is rather 
simple. The electronically active valence band consists 
predominantly of 5p orbitals of the Te atoms from the Te 
planes. The only significant role in the formation of the 
CDW is played by the perpendicular chains of in-plane 
bpx and 5py orbitals, since the energy of completely filled 
bpz is pushed below the Fermi level as indicated by the 
first principle band structure calculations 0, [l^ . 



II. STM MEASUREMENTS 

We have performed low-temperature STM topography 
and spectroscopy of CeTes at 77 K to characterize the 
CDW state. CeTes crystals were grown by a halide flux 
method, as described in Ref. [17]. The crystals were care- 
fully cleaved with adhesive tape and quickly placed in a 
vacuum system for subsequent STM measurements at a 
temperature of 77 K. A representative unfiltered STM 
image is shown in Fig. [2l Both the net of Te atoms and 
the CDW modulation are clearly visible, as indicated. 
The CDW modulation is oriented at 45° to the Te net. 

Fig. [3] (a) shows the Fourier transform of the 77 K 
real-space data obtained from 16 24x24 nm images us- 
ing a straightforward averaging procedure The Te 
square lattice peaks are labeled L. Peaks related to the 
CDW are oriented 45° clockwise to the square lattice 



peaks. To examine the peaks carefully, we take a line cut 
along the CDW direction in Fig. [3] (b). We expect the 
CDW peak to be located near qcDW ^ 2/7 x 27r/c = 
4.1 nm~^. Indeed, we find a prominent peak at 3.9 nm~^, 
which we label qcDW in bold. The peak at 14.3 nm~^, 
labeled qatom^ also appears in the direction perpendicular 
to the CDW. Hence it is consistent with a larger square 
lattice oriented at 45° to the surface Te net. As seen 
in Fig. [H this pattern is consistent with either the first 
subsurface Ce layer or the first subsurface Te layer. A 
similar Fourier transform peak was observed by Fang et 
al.^ while working with the related material TbTes, who 
also attributed it to sensitivity to the subsurface struc- 
ture [i^. Surprisingly, qatom and qcDW have roughly 
the same magnitude, which is approximately equal to 
the magnitude of the Te-net Fourier peaks, labeled L in 
Fig. [3] (a). Theoretical calculations given in section [Vl 
will address the relative magnitudes of the contributions 
of the tunneling current from the surface Te net and sub- 
surface Ce and Te atoms. 

Because the Fourier transform lacks phase informa- 
tion, the square lattice corresponding to peak qatom could 
possibly be attributed to the apparent dimerizations of 
the surface Te net reported by Fang and coworkers while 
working with the related material TbTea [10]. However, 
we observe no evidence of dimers in topographic images 
of CeTea when reproducing the same tunneling condi- 
tions (Fig. 2]). Moreover, we have not seen evidence for 
dimers in the hundreds of images we have taken of CeTes 
over a broad set of bias voltages in the range of ±800 mV. 
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FIG. 3: (a) The Fourier transform of the low temperature 
STM data. Enhanced noise along the vertical axis is an ar- 
tifact due to the scan direction. Horizontal and vertical axes 
are wave vector components kx and ky. The square Te net 
gives rise to four distinct peaks (L). Peaks at 45° to Te net 
are consistent with the CDW peaks. The CDW peak qcnw, 
as well as peak qatom, which is consistent with underlying 
atomic structure, are labeled, (b) Subset of the data from 
the Fourier transform along a path from the origin in the di- 
rection of the CDW. Noise in the Fourier transform becomes 
significantly larger near the origin. This is due to the impu- 
rities present in the real-space data. As the guide to the eye, 
the dashed line indicates the background noise. 



III. WAVE VECTOR MIXING 

Returning to Fig. [3] (b), four peaks are present in ad- 
dition to QcDW and qatom- In our original study, we 
interpreted similar structure as satellite peaks, which 
supported the interpretation of a discommensurated 
CDW [9[. However, Fisher has suggested an alternative 
wave- vector- mixing explanation, which we explore here 
in detail [11]. 

Consider an STM signal acquired on a surface with 
a uniformly incommensurate CDW and lattice modula- 
tion, both of which are described by sinusoidal waves 
with wave vectors kcDW and katom along the x-direction. 
The Fourier transform along this direction would ide- 
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FIG. 4: STM topography of CeTes at a sample bias of - 
800 mV and tunneling current of -0.05 nA, at room tempera- 
ture in ambient conditions, the same tunneling conditions as 
Fig. 2(a) of reference [13]. We do not observe evidence of 
dimerizations, in contrast to the work of Fang et ai with the 
related material, TbTes. These data were filtered to reduce 
noise. 



ally exhibit exactly two peaks, one at kcDW and one 
at katom- Implicit in the ideal case is the assumption 
that the STM signal is proportional to the sum of the 
two sine waves. However, in reality the signal may have 
a significant component that resembles the product of 
the two waves. In reference [18], Tomic showed explicitly 
that this is indeed the case if an asymmetry existed in 
the way the CDW couples to the peaks of the atomic 
signal compared to the troughs. In other words, the 
signal will exhibit a contribution similar to the prod- 
uct of the two waves if the CDW signal is stronger at 
the atomic lattice sites and weaker at the locations be- 
tween atoms, or vice versa. Algebraically, the product 
of two sine waves can be expressed as a sum and a 
difference: sin{kcDwx)sin{katomX) = \cos{kcDwx — 

katomX) - \cOs{kcDWX + katomX) . HcUCC WC CXpCCt 

this effect to give extra peaks in the Fourier transform 
at kcDW + katom and kcDW - katom- More generally, 
Tomic showed that such a peak-trough asymmetry leads 
to wave- vector mixing, additional Fourier peaks at linear 
combinations of the two wave vectors. 

With this in mind, we interpret the additional peaks 
in the FT of our STM data as linear combinations of 
QCDW and qatom- In Fig-[3](b), the additional peaks are 
labeled as qatom - '^Qcdw, '^Qcdw, Qatom - qcDW and 
Qatom + QCDW- The arrows mark the precise wave vector 
value at which we would expect each mixed peak to occur; 
indeed we see excellent agreement between the arrows 
and the actual peak locations. This analysis is similar 
to the analysis presented by Fang et al to describe the 
TbTes data (lO|. 

As all of the clear Fourier peaks are accounted for 
without invoking discommensurations, our STM data of 
CeTea are consistent with a uniformly incommensurate 
CDW. However, this analysis does not preclude the pres- 
ence of discommensurations. It is possible that some 
of the unlabeled features such as the small peaks near 
16.5 nm~^ and 20 nm~^ are satellite peaks indicative of 
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FIG. 5: Local DOS around the Fermi level for CeTes. The 
Fermi level corresponds to zero bias voltage. The estimated 
CDW gap size in CeTes is about 360 meV, as indicated by ver- 
tical arrows. The spectroscopy was carried out when the tip 
was located directly above a Te atom. The effect of thermal 
smearing of the data is approximately 3.5 ksT ^ 25 meV at 
77 K temperature, hence it is not significant over the plotted 
range. 



discommensurations [9]. However these features are just 
above the noise level of the data. 

Given that CDW and lattice modulations are typically 
well-described by undistorted sinusoidal waves, it is sur- 
prising that the wave- vector- mixing effect is so large. In 
light of the analysis by Tomic and our assertion that the 
lattice signal in this case corresponds to the subsurface 
Ce or Te, we conjecture the following. The CDW modu- 
lates the amplitude of the tunneling signal of the surface 
Te atoms; these are located at trough positions with re- 
spect to the subsurface lattice. But the CDW is expected 
to only weakly couple to the subsurface atoms them- 
selves. We believe that this is the source of the peak- 
trough asymmetry. Unfortunately, the effect is rather 
subtle with respect to the real-space patterns; given the 
noise in the measurement, it is difficult to confirm this 
conjecture by examining the direct images. 



IV. SPECTROSCOPY 

Formation of the CDW state is expected to be ac- 
companied by an energy gap opening up at the Fermi 
level. Using the point spectroscopy mode we can probe 
the CDW gap at different locations of the sample surface 
and estimate its size. Spectra were acquired at a tem- 
perature of 77 K at various locations of the tip above the 
exposed Te plane on the surface of CeTes. 

Fig. [5] shows the characteristic density of states as ob- 
tained when the tip is above a Te atom. The data repre- 
sent an average of 168 measurements performed consecu- 
tively. Arrows mark the approximate edges of the CDW 



gap, estimated to have a width of 360 meV. The mag- 
nitude and shape of the dl/dV curves are remarkably 
similar to the previous TbTes measurement by Fang and 
coworkers and may be compared to the ARPES results of 
400 meV [7, iO] . The local density of states has a shape 
that is suggestive of subgap states, characterized by the 
non-zero V-shaped structure within the gap. 



V. THEORETICAL STUDIES 

A. Method of Calculation 

The STM data obtained for CeTes was simulated us- 
ing electronic structure calculations within density func- 
tional theory (DFT) 0. It is well known that the 
local (spin) density approximation L(S)DA fails to de- 
scribe the correct ground state of a systems contain- 
ing transition metal or rare-earth metal atoms. L(S)DA 
always puts the partially filled d or f bands right at 
the Fermi level {Ep)^ predicting metallic character with 
itinerant d or f electrons, which is obviously not cor- 
rect. Strong Coulomb repulsion between localized d 
(or /) electrons suppresses the charge fluctuations in- 
herent in a metallic system. In order to describe cor- 
rectly the ground state of such systems, one has to go 
beyond standard L(S)DA, and take into account the 
strong electron-electron correlations. One of the success- 
ful approaches is the L(S)DA+U method |2l|, |22|, l23|, l2| , 
in which the localized d or f electrons and the delo- 
calized s and p electrons are treated differently. The 
orbital-dependent Coulomb potential is only taken into 
account for localized states, while the delocalized states 
are treated by orbital-independent L(S)DA type poten- 
tial. For CeTes the LSDA+U approach was used based 
on full potential linearized augmented plane wave + lo- 
cal orbital (FPLAPW+lo) method [ll, |26| as imple- 
mented in the Wien2k package [27^]. The on-site elec- 
tron correlation was taken into account for the Ce /- 
states with Ueff{Ce) = 6.8 eV. For the exchange and 
correlation functional the local spin density approxima- 
tion (LSDA) [2^ was used. The value of the convergence 
parameter RKmaxi which is defined as the product of the 
minimal atomic sphere radius {R) and the largest recip- 
rocal lattice vector {Kmax) of the plane wave basis, was 
chosen as RKmax = 7. We use 2.5 a.u. for the muffin-tin 
radii of all Ce and Te atoms. The Brillouin zone (BZ) 
was sampled by a dense mesh of 625 irreducible k-points 
in the /c^ = plane. Spin-orbit (SO) interaction was 
included using the second variational treatment [2^, [sOl • 

B. STM Simulation 

The CeTes surface was modeled by a periodic slab 
geometry separated by a vacuum region of 14 A. The 
slabs consisted of repeating supercells, each containing 
two unit cells along the crystallographic long axis (b-axis 
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FIG. 6: To simulate the STM image, we give the contour plot 
of the charge density projected onto a plane at 3.0 A above 
the surface, in the energy range of 0.0 — 0.1 eV. The brighter 
(white) spots come from the surface Te. The subsurface Ce 
atoms contribute more to the charge density than the sub- 
surface Te atoms. The Fermi level (^f) corresponds to 0.0 
eV. Moreover, the subsurface modulation, considering only 
the regions between surface Te, is of comparable magnitude 
to the signal from the surface Te net. 



for CeTes). The supercell was constructed using the the- 
oretical lattice constants of the bulk CeTes. The thick- 
ness of the slab constructed this way is 47.5 A, which 
is sufficient to eliminate the spurious surface-surface in- 
teraction. Geometry optimization was performed on the 
surface Te square lattice and on the Ce and Te layers 
closest to the surface. The calculations were carried out 
on the high symmetry structure. Since calculations us- 
ing supercell models for the incommensurate CDW are 
not feasible, we focus on the nature of the electronic wave 
functions near Ef^ on the undistorted structure, with the 
goal of elucidating the contributions from the subsurface 
atoms to these wave functions. 

According to the theory of tunneling between a real 
solid surface and model probe tip, the tunneling current 
in the first-order perturbation theory is given by [31, 32} : 



/ = 



27re 



(1) 

where f{E) is the Fermi function, V is the applied volt- 
age, and Ej^ are the energies of the states ip^ and ipj^^ 
of the probe tip and the surface, respectively. M^^ is the 
tunneling matrix element between ip^ and tpj^. Taking 
the limit of small voltage (100 mV) and low temperature 
(77 K) and approximating the tip with a spherical wave 
function, the tunneling current can be written as: 
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FIG. 7: The only significant contribution to the DOS near Ep 
comes from the surface Te px and py states, (a) Partial DOS 
associated with the px and py orbit als of the surface Te atoms. 
The flat DOS near is characteristic to the CDW systems, 
(b) The pz orbital of the surface Te is occupied and located 
mostly below Ef- (c),(d) The p-states of the subsurface Te 
atoms are occupied, because the trivalent Ce provides two 
electrons within the Ce-Te slab. 



where Ef is the Fermi energy. When Ei, = Ef^ the 
tunneling current is proportional to the local density of 
states (LDOS) at the position of the tip (ro), as given by 
equation O Using the slab geometry described above we 
simulate a constant-height-mode STM image by calculat- 
ing the charge density in an energy range of 0.0 — 0.1 eV 
around Ef taken at a distance of 3.0 A above the surface 
Te net, where 0.0 eV is the Fermi energy. Fig. [6] shows 
the contour plot of this charge density. The main con- 
tribution comes from the surface Te atoms arranged in a 
square lattice. The subsurface Ce and Te atoms are ar- 
ranged in puckered double layers in which the Ce and Te 
atoms form square lattices oriented at 45° to the surface 
Te net. As shown in Fig.[6l the subsurface Ce contributes 
more to the charge density than the subsurface Te, sug- 
gesting that tunneling is more likely to appear from the 
Ce sites than from the subsurface Te sites. This is con- 
sistent with the experimental finding that the peak q at 
14.3 nm~^ in Fig. [3] corresponds to the underlying struc- 
ture closest to the surface, namely to the square lattice 
of Ce atoms. 

We have obtained similar results when we calculated 
the LDOS in an energy range of ±0.08 eV, at distances 
of 2.5 A and 3.0 A above the surface Te net. 



C. Electronic Structure 

Figure [71 shows the calculated partial density of states 
(PDOS) associated with the p-orbitals of the surface and 
subsurface Te atoms. Since the spin-up and spin-down 
DOS associated with the Te atoms are identical, we only 
show the spin-up channel. As pointed out by DiMasi et 
al. [H, [s^ in RETea compounds the rare earth is triva- 
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FIG. 8: Spin-polarized, partial DOS associated with the sub- 
surface Ce atom. The peak at -2.5 eV in the spin-up DOS 
corresponds to one occupied f- level while the peaks between 
3.0-5.0 eV (spin-up and spin-down) represent the 13 unoc- 
cupied f-states. The inset shows that the DOS near Ep is 
nonzero. 



lent, it provides two electrons for bonding to the puck- 
ered double-layer of RE-Te and one electron to the Te 
square net. Therefore we expect the p-states of the sub- 
surface Te atoms (which are part of the RE-Te slabs) 
to be completely occupied, while the p-states surface Te 
atoms (which form the square net) should be partially 
filled. As shown in Fig. [7J our electronic structure calcu- 
lation agrees with this observation and with the earlier 
theoretical results [l5|, ll6|: the only significant contri- 
bution to the DOS at Ep comes from the Px and Py 
orbitals of the surface Te. The partial DOS associated 
with these two orbitals (Fig. [71 (a)) are rather flat at 
Ef, which is a characteristic of the CDW systems in the 
high symmetry phase. The Pz orbital of the surface Te 
atom (Fig. [3 (b)) and the the p-states of the subsurface 
Te atoms (Fig. [7](c), (d)) are nearly fully occupied and 
located mostly below Ep- 

In order to understand why the subsurface Ce atoms 
contribute so much to the STM charge density calcu- 
lation, we give the spin polarized DOS associated with 
the subsurface Ce atom in Fig. [HI The spin-up states 
are represented in the positive (upper) region, while the 
spin-down states are shown in the negative (lower) re- 
gion of the graph. The main contribution to the DOS 
comes from the Ce /-states: the narrow, sharp peak in 
the spin-up DOS located at ~ 2.5 eV below Ep represents 
the only one occupied Ce /-level. The broader peaks lo- 
cated between 3.0 eV and 5.0 eV above Ef^ in both the 
spin- up and spin-down DOS, come from the rest of the 
13 empty Ce /-states (6 spin- up and 7 spin-down states). 
The splitting between the occupied and empty states is 
approximately equal to the chosen value of the Coulomb 
repulsion within the /-shell Ueff{Ce) = 6.8 eV. The inset 
of Fig. [8] gives the DOS associated with the Ce orbitals 
near Ep. There is a small but finite contribution coming 



from both the hybridization of /-states with Te p-bands 
and Ce d states. These states will contribute to tunneling 
measurements above the Ce sites. 

Now let us discuss the tunneling results in the light 
of our theoretical calculations. As discussed in sec. IV B[ 
the STM tunneling current is proportional to the LDOS 
at the position of the tip. In reality however there is a 
matrix element between the states from which the tun- 
neling occurs (host) and the tip state. Assuming that 
the tip state can be approximated by a smooth 5-like 
function, the symmetry of the host state will determine 
the strength of the tunneling current. For example, for 
the surface Te, Px and Py states that make dominant 
contribution to the DOS near Ep will have zero (very 
small) contribution to the tunneling current. As dis- 
cussed above, the intensity of the FT of the STM signal 
as a function of q along the direction of the CDW wave- 
vector (which is rotated by 45° with respect to the direc- 
tion of the square lattice peaks) gives peaks at qcDW = 
3.9 nm~^ and at Qatom = 14.3 nm~^, the latter is ascribed 
to the subsurface Ce or Te lattice. From our PDOS cal- 
culations we find that for the subsurface Te, the Px and 
Py orbitals contribute ~ 0.004 stats/eV/spin near EF, 
whereas the Pz orbital contributes ~ 0.02 states/eV/spin. 
Thus both due to symmetry and small PDOS, the Px and 
Py states will not contribute to the tunneling current, 
only the Pz orbital will contribute. The PDOS associated 
with the Ce orbitals are ~ 0.03 states/eV/spin compa- 
rable to the subsurface Te Pz . However the lattice of Ce 
atoms is a distance of 2.53 A from the surface Te net, 
whereas the subsurface Te lattice is a distance of 3.44 A. 
Hence the Ce atoms are nearly 1 A closer to the STM 
tip; it is therefore likely that the tunneling current above 
Ce sites is larger than that above the subsurface Te sites. 
Furthermore, since the surface Te Pz contribution to the 
DOS near Ej? is rather small (~ 0.004 stats/eV/spin) 
our charge density analysis suggests that the tunneling 
currents above the Ce sites and the surface Te net sites 
are comparable. This is the origin of the surprisingly high 
sensitivity to the subsurface Ce, consistent with both the 
experiment, (Fig. [3]) and theoretical simulation (Fig. [6]). 
This result is significant because it may be possible to 
probe the nature of Ce / states through careful tunnel- 
ing measurements. 



VI. SUMMARY 

STM constant-current-mode images and spectroscopy 
measurements of CeTea were acquired at 77 K. The ef- 
fects of the CDW were clearly resolved in both the images 
and the spectroscopy curves. In addition to the CDW, 
the images show the expected signal from the surface 
net of Te atoms and a large contribution from subsurface 
structure of approximately the same magnitude. Con- 
trary to a study by Fang et al. [10], no evidence of 
dimerization of the surface Te net was observed. 

Fourier transform analysis of the STM images showed 
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two principal peaks along the direction of the charge den- 
sity wave, one from the CDW at Qcdw =3.9 nm~^ and 
the other from a subsurface lattice at qatom = 14.3 nm~^. 
In addition to these peaks four others were observed 
which we show are well described by a wave-vector- 
mixing effect of the principle wave vectors. We believe 
the high degree of wave vector mixing is caused essen- 
tially by the fact that the CDW, which exists in the sur- 
face layer, distorts preferentially the troughs between the 
subsurface atoms. With regard to the uniformity of the 
CDW, our present study does not provide support for 
discommensurations, in contrast to the interpretation of 
our earlier STM work [9]. However, the data presented 
here do not rule out discommensurations and the issue 
remains an open question. 

To better understand the nature of the tunneling sig- 
nal and the large contribution from subsurface atoms, we 
theoretically considered the symmetry of the tunneling 
matrix elements and performed density function theory 



calculations. We concluded that the dominant contribu- 
tion to the subsurface signal is from a lattice of Ce atoms 
2.53 A below the surface Te net. Moreover a simulated 
STM image constructed from our calculations confirms 
that the modulation of the signal arising from this lat- 
tice is comparable in magnitude to the signal from the 
surface Te net. 
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